QCD thermodynamics and confinement from 
a dynamical quasiparticle point of view 



W. Cassing a '* 



^Institut fur Theoretische Physik, Universitat Giessen, Heinrich-Buff-Ring 16, 

D-35392 Giessen, Germany 



Abstract 

In this study it is demonstrated that a simple picture of the QCD gluon liquid 
emerges in the dynamical quasiparticle model that specifies the active degrees of 
freedom in the time-like sector and yields a potential energy density in the space-like 
sector. By using the time-like gluon density (or scalar gluon density) as an indepen- 
dent degree of freedom - instead of the temperature T as a Lagrange parameter - 
variations of the potential energy density lead to effective mean-fields for time-like 
gluons and an effective gluon-gluon interaction strength at low density. The latter 
yields a simple dynamical picture for the gluon fusion to color neutral glueballs when 
approaching the phase boundary from a temperature higher than T c and paves the 
way for an off-shell transport theoretical description of the parton dynamics. 
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1 Introduction 



The formation of a quark-gluon plasma (QGP) and its transition to interacting 
hadronic matter - as occurred in the early universe - has motivated a large commu- 
nity for several decades (cf. [1] and Refs. therein). Early concepts of the QGP were 
guided by the idea of a weakly interacting system of partons (quarks, antiquarks 
and gluons) since the entropy s and energy density e were found in lattice QCD to 
be close to the Stefan Boltzmann (SB) limit for a relativistic noninteracting system 
[2]. However, this notion had to be given up in the last years since experimental 
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observations at the Relativistic Heavy Ion Collider (RHIC) indicated that the new 
medium created in ultrarelativistic Au+Au collisions was interacting more strongly 
than hadronic matter. Moreover, in line with earlier theoretical studies in Refs. [3- 
5] the medium showed phenomena of an almost perfect liquid of partons [6,7] as 
extracted from the strong radial expansion and elliptic flow of hadrons as well the 
scaling of the elliptic flow with parton number etc. The latter collective observables 
have been severely underestimated in conventional string/hadron transport models 
[8-10], but hydrodynamical approaches did quite well in describing (at midrapidity) 
the collective properties of the medium generated during the early times for low and 
moderate transverse momenta [11,12]. Soon the question came up about the con- 
stituents of this liquid; it might be some kind of i) "epoxy" [13], i.e. a system of 
resonant or bound gluonic states with large scattering length, ii) a system of chirally 
restored mesons, instanton molecules or equivalently giant collective modes [14], iii) 
a system of colored bound states of quarks q and gluons g, i.e. gq, qq, gg etc. [15], iv) 
some 'string spaghetti' or 'pasta' etc. In short, many properties of the new phase are 
still under debate and practically no dynamical concepts are available to describe 
the freezeout of partons to color neutral hadrons that are subject to experimental 
detection. 

Lattice QCD (1QCD) calculations provide some guidance to the thermodynamic 
properties of the partonic medium close to the transition at a critical temperature 
T c up to a few times T c , but 1QCD calculations for transport coefficients presently 
are not accurate enough [16] to allow for firm conclusions. Furthermore, it is not 
clear whether the partonic system really reaches thermal and chemical equilibrium 
in ultrarelativistic nucleus-nucleus collisions and nonequilibrium models are needed 
to trace the entire collision history. The available string/hadron transport models 
[17-19] are not accurate enough - as pointed out above - nor do partonic cascade 
simulations [20-23] (propagating massless partons) sufficiently describe the reaction 
dynamics when employing cross sections from perturbative QCD (pQCD). This also 
holds - to some extent - for the Multiphase Transport Model AMPT [24] since it 
includes only on-shell massless partons in the partonic phase as in Ref. [21]. The 
same problem comes about in the parton cascade model of Xu and Greiner [25] where 
additional 2<-^ 3 processes like gg <-> ggg are incorporated. On the other hand it is 
well known that strongly interacting quantum systems require descriptions in terms 
of propagators D with sizeable selfenergies II for the relevant degrees of freedom. 
Whereas the real part of the selfenergy gives contributions to the energy density, the 
imaginary parts of II provide information about the lifetime and/or reaction rate 
of time-like 'particles' [4]. In principle, off-shell transport equations are available in 
the literature [26-28], but have been applied only to dynamical problems where the 
width of the quasiparticles stays moderate with respect to the pole mass [29]. On 
the other hand, the studies of Peshier [30,31] indicate that the effective degrees of 
freedom in a partonic phase should have a width 7 in the order of the pole mass M 
already slightly above T c . 

The present study addresses essentially three questions: i) Do we understand the 
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QCD thermodynamics in terms of dynamical quasiparticles down to the phase bound- 
ary in a 'top down' scenario and what are the effective degrees of freedom as well 
as energy contributions? ii) Can such a quasiparticle approach help in defining an 
off-shell transport model that - at least in thermal equilibrium - reproduces the 
thermodynamic results from 1QCD? iii) Are there any perspectives in modeling the 
transition from partonic to hadronic degrees of freedom in a dynamical way? 

The present work is exploratory in the sense that it is restricted to a pure gluonic 
system of N% — 1 gluons with two transverse polarisations, i.e. degeneracy d g = 16 for 
the gluonic quasiparticles that are treated as relativistic scalar fields. Note, however, 
that the qualitative features stay the same when adding light quark degrees of free- 
dom [31]; this finding is well in line with the approximate scaling of thermodynamic 
quantities from 1QCD when dividing by the number of degrees of freedom and scaling 
by the individual critical temperature T c which is a function of the different number 
of parton species [32]. 

The outline of the paper is as follows: After a short recapitulation of the dynamical 
quasiparticle model in Section 2 new results on the space-like and time-like parts 
of observables are presented that allow for a transparent physical interpretation. In 
Section 3 we will examine derivatives of the space-like part of the quasiparticle energy 
density with respect to the time-like (or scalar) density which provides information on 
gluonic mean fields and their effective interaction strength. The implications of these 
findings with respect to an off-shell transport description are pointed out throughout 
the study. A summary and extended discussion closes this work in Section 4. 



2 Off-shell elements in the DQPM 

2. 1 Reminder of the DQPM 

The Dynamical QuasiParticle Model (DQPM) 1 adopted here goes back to Peshier 
[30,31] and starts with the entropy density s in the quasiparticle limit [33], 

where n{uj/T) = (exp(cj/T) — 1) _1 denotes the Bose distribution function, A stands 
for the scalar quasiparticle propagator and II for the quasiparticle selfenergy which is 
considered here to be a Lorentz scalar. In principle, the latter quantities are Lorentz 
tensors and should be evaluated in a nonperturbative framework. However, a more 
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practical procedure is to use a physically motivated Ansatz with a Lorentzian spectral 
function, 




P{UJ) ~ - I ( W -£)2 +7 2 (c , + £)2 +72 )> ( 2 ) 



and to fit the few parameters to results from 1QCD. With the convention E 2 (p) = 
p 2 + M 2 — 7 2 , the parameters M 2 and 7 are directly related to the real and imaginary 
parts of the corresponding (retarded) self-energy, II = M 2 — 2i>yuj. It should be 
stressed that the entropy density functional (1) is not restricted to quasiparticles of 
low width 7 and thus weakly interacting particles. In fact, in the following it will 
be shown that a novel picture of the hot gluon liquid emerges because 7 becomes 
comparable to the quasiparticle mass already slightly above T c [30,31]. 

Following [34] the quasiparticle mass (squared) is written in (momentum-independent) 
perturbative form, 

M 2 (T) = fg 2 T 2 , (3) 



with a running coupling (squared), 

g 2 (T/T c ) = niVcln ( A 2( T / Tc _ Ts / Tc ) 2 ' ( 4 ) 



which permits for an enhancement near T c [34,35]. It will be shown below that an 
infrared enhancement of the coupling - as also found in the 1QCD calculations in 
Ref. [36] for the long range part of the q — q potential - is directly linked to the 
gluon fusion/clustering scenario. In order to quantify this statement the coupling 
a s (T) = g 2 (T) / (Att) is shown in Fig. 1 as a function of T/T c in comparison to the long 
range part of the strong coupling as extracted from Ref. [36] from the free energy of a 
quark- ant iquark pair in quenched 1QCD. For this comparison the actual parameters 
A = 2.42, T s /T c = 0.46 have been adopted as in Ref. [4]. The parametrization (4) is 
seen to follow the 1QCD results - also indicating a strong enhancement close to T c - 
as a function of temperature reasonably well. One should recall that any extraction 
of coupling constants a s (T) from 1QCD is model dependent and deviations from (or 
agreement with) lattice 'data' have to be considered with care. The argument here is 
that the specific 'parametric form' of Eq. (4) is not in conflict with 1QCD and that 
the coupling a s and consequently the quasiparticle mass M(T) has the right order 
of magnitude. 



4 




0.0 i i i i i i i i i 

1 2 3 456789 10 



T/T c 

Fig. 1 . The coupling a s (T) = g 2 (T) / (47r) (solid red line) as a function of T/T c in comparison 
to the long range part of the strong coupling as extracted from Ref. [36] from the free energy 
of a quark-antiquark pair in quenched 1QCD (for N T = 8). 

The width 7 is adopted in the form 7 ~ g 2 Thig~ l [37] or, equivalently, in terms of 
M [30], as 

3 MHJ^ c 
7(T) = 4^^^ Tln (M(7W (5) 



where c = 14.4 (from [4]) is related to a magnetic cut-off. In case of the pure Yang- 
Mills sector of QCD the physical processes contributing to the width 7 are both 
99 99 scattering as well as splitting and fusion reactions gg «-> g or gg <-> ggg, 
999 9999 e ^ c - Note that the ratio j(T)/M(T) ~ ^ln(c/^ 2 ) approaches zero only 
asymptotically for T — > 00 such that the width of the quasiparticles is comparable to 
the mass for all practical energy scales on earth; the ratio j(T)/M(T) drops below 
0.5 only for temperatures T > 1.25 • 10 5 T c (for the parameters given above). 

For the choice (2) for the spectral function the scalar effective propagator reads, 



which can easily be separated into real and imaginary parts. The entropy density (1) 
then reads explicitly [31], 

s dqp (T) = d 9 j^{~ HI - e^ T ) + ^n{u p /T)) 

r du d 3 p dn ( 2 7 ^ - u 2 ) \ 

+ 9 J 27r(27r)3 dT\ aC [ lu 2 -lu 2) (^-lu^ + A^lu 2 ) ' [<) 
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using uj p = Vp 2 + M' 2 . The first line in (7) corresponds to the familiar on-shell quasi- 
particle contribution s while the second line in (7) corresponds to the contribution 
originating from the finite width 7 of the quasiparticles and is positive throughout 
but subleading (see below). 

The pressure P now can be evaluated from 
OP 



by integration of s over T, where from now on we identify the 'full' entropy density s 
with the quasiparticle entropy density s dqp . Note that for T <T C the entropy density 
drops to zero (with decreasing T) due to the high quasiparticle mass and the width 7 
vanishes as well because the interaction rate in the very dilute quasiparticle system 
becomes negligible. Since the pressure for infinitely heavy (noninteracting) particles 
also vanishes the integration constant for the pressure P - when integrating (8) - 
may safely be assumed to be zero, too. 

The energy density e then follows from the thermodynamical relation [34,38] 

e = Ts-P (9) 

and thus is also fixed by the entropy s(T) as well as the interaction measure 

W{T) := e(T) - 3P(T) =Ts-4P (10) 



that vanishes for massless and noninteracting degrees of freedom. 

In Ref. [4] a detailed comparison has been presented with the lattice results from 
Ref. [39] for the pure gluonic sector to the quasiparticle entropy density (7) for the 
parameters given above. The agreement with the lattice data is practically perfect 
[4,30]. Needless to point out that also P(T),e(T) and W(T) well match the lattice 
QCD results for 1 < T/T c < 4 [4,31] due to thermodynamical consistency. The same 
parameters are also adopted for the following calculations. 



2.2 Time-like and space-like quantities 



For the further argumentation it is useful to introduce the shorthand notation 

Trp • • • = d 9 J^^- 3 2u p{u) 6H n{u/T) 0(±P 2 ) • • • (11) 

with P 2 = u) 2 — p 2 denoting the invariant mass squared. The 0(±P 2 ) function in 
(11) separates time-like quantities from space-like quantities and can be inserted for 
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any observable of interest. 



As the first quantity we consider the entropy density (7). Its time-like contribution 
is almost completely dominated by the first line in (7) - that corresponds to the on- 
shell quasiparticle contribution s - but also includes a small contribution from the 
second line in (7) which is positive for T below about 1.5 T c and becomes negative for 
larger temperature. This time-like part s + is shown in Fig. 2 by the dotted blue line 
(multiplied by (T c /T) 3 ). The second line in (7) - as mentioned above - corresponds 
to the contribution originating from the finite width 7 of the quasiparticles and also 
has a space-like part s~ which is dominant (for the second line in (7)) and displayed 
in Fig. 2 by the lower red line (multiplied by (T c /T) 3 ). Though s~ is sub leading in 
the total entropy density s = s + + s~ (thick solid green line in Fig. 2) it is essential 
for a proper reproduction of s(T) close to T c (cf. [31]). Note that the total entropy 
density s is not very different from the Stefan Boltzmann entropy density ssb for 
T > 2T C as shown in Fig. 2 by the upper thin line (multiplied by (T c /T) 3 ). 
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Fig. 2. The time-like contribution to the entropy density s + (dotted blue line), the space-like 
contribution s~ (lower red line) and the total entropy density s = s + + s~ (thick solid green 
line) as a function of T/T c . All quantities have been multiplied by the dimensionless factor 
(T c /T) 3 ) assuming T c = 0.26 GeV for the pure gluonic system [40]. The upper solid black 
line displays the Stefan Boltzmann limit ssb for reference. 



Further quantities of interest are the quasiparticle 'densities' 
A^ ± (T) 



Tr* 1 



(12) 



that correspond to the time-like (+) and space-like (-) parts of the integrated distri- 
bution function. Note that only the integral of iV + over space has a particle number 
interpretation. In QED this corresponds to time-like photons (7*) which are virtuell 



7 



in intermediate processes but can also be seen asymptotically by dileptons (e.g. e + e 
pairs) due to the decay 7* — > e + e~ [17]. 



A scalar density N s , which is only defined in the time-like sector, is given by 

( \fP 1 \ 
N S (T) = Tr+ V ■ 



UJ 



(13) 



and has the virtue of being Lorentz invariant. Moreover, a scalar density can easily 
be computed in transport approaches for bosons and fermions [17,41] which is of 
relevance for the argumentation in Section 3. 
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Fig. 3. Upper part: The scalar density N s (lower orange line), the time-like density N + 
(blue line), the space- like quantity N~ (red line) and the sum N = N + + N~ (thick solid 
green line) as a function of T/T c assuming T c = 0.26 GeV for the pure gluonic system 
[40]. The upper solid black line displays the Stefan Boltzmann limit N$b for reference. All 
quantities are multiplied by the dimensionless factor (T c /T) 3 . Lower part: The ratio of the 
scalar density N s to the time-like density iV + as a function of the scaled temperature T/T c . 



The actual results for the different 'densities' (multiplied by (T c /T) 3 ) are displayed 
in the upper part of Fig. 3 where the lower orange line represents the scalar density 
N s , the blue line the time-like density iV + , the red line the space-like quantity N~ 
and the thick solid green line the sum N = N + + N~ as a function of T/T c assuming 
(as before) T c = 0.26 GeV for the pure gluonic system [40]. It is seen that N + is 
substantially smaller than N~ in the whole temperature range up to 10 T c where it 
is tacitly assumed that the DQPM also represents 1QCD results for T > 4T C , which 
is not proven explicitly, but might be expected due to the proper weak coupling 
limit of (3), (5) (cf. Fig. 1). The application of the DQPM to 10 T c is presented 
in Fig. 3 since the initial state at Large Hadron Collider (LHC) energies might be 
characterized by a temperature above 4 T c ; note that the properties of the partonic 
phase will be explored from the experimental side in the near future at LHC. Quite 
remarkably the quantity iV follows closely the Stefan Boltzmann limit N$b for a 
massless noninteracting system which is given in Fig. 3 by the upper thin solid line 
and has the physical interpretation of a gluon density. Though N differs by less than 
15% from the Stefan Boltzmann (SB) limit for T > 2Tc the physical interpretation 
is essentially different! Whereas in the SB limit all gluons move on the light cone 
without interactions only a small fraction of gluons can be attributed to quasiparticles 
with density N + within the DQPM that propagate within the lightcone. The space- 
like part N~ corresponds to 'gluons' exchanged in t-channel scattering processes 
and thus cannot be propagated explicitly in off-shell transport approaches without 
violating causality and/or Lorentz invariance. 

The scalar density N s follows smoothly the time-like density N + as a function of 
temperature which can be explicitly seen in the lower part of Fig. 3 where the ratio 
N s /N + is shown versus T/T c . Consequently, the scalar density N s uniquely relates 
to the time-like density N + or the temperature T in thermal equilibrium which will 
provide some perspectives for a transport theoretical treatment (see Section 3). 

The separation of N + and N~ so far has no direct dynamical implications except 
for the fact that only the fraction N + can explicitly be propagated in transport as 
argued above. Thus we consider the energy densities, 



that specify time-like and space-like contributions to the quasiparticle energy density. 
It is worth pointing out that the quantity T 00 = Too + ^oo i n case °f a conventional 
quasi-particle model with vanishing width 7 in general is quite different from e in (9) 
because the interaction energy density in this case is not included in (14), i.e. 



Too = T+ = d g du—^ 2cu 5(u 2 — M 2 — p 2 ) G(u) 0(±P 2 ) n(u/T) to (15) 



7*(T) = tr ± u , 



(14) 




since uj 2 — p 2 = 



M 2 = P 2 > due to the mass-shell 5-function. 
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How does the situation look like in case of dynamical quasiparticles of finite width? 
To this aim we consider the integrand in the energy density (14) which reads as (in 
spherical momentum coordinates with angular degrees of freedom integrated out) 

I(uj,p) = ip 2 uj 2 p(uj,p 2 )n(uj/T). (16) 



Here the integration is to be taken over uj and p from to oo. The integrand I(u, p) is 
shown in Fig. 4 for T = 1.02T C (l.h.s.) and T = 2T C (r.h.s.) in terms of contour lines. 
For the lower temperature the gluon mass is about 0.91 GeV and the width 7 « 0.15 
GeV such that the quasiparticle properties are close to a p-meson in free space. In 
this case the integrand I(u,p) is essentially located in the time-like sector and the 
integral over the space-like sector is subdominant. This situation changes for T = 2T C 
where the mass is about 0.86 GeV while the width increases to 7 ~ 0.56 GeV. As one 
observes from the r.h.s. of Fig. 4 the maximum of the integrand is shifted towards 
the line uo = p and higher momentum due to the increase in temperature by about 
a factor of two; furthermore, the distribution reaches far out in the space-like sector 
due to the Bose factor n(u/T) which favors small to. Thus the relative importance of 
the time-like (+) part to the space-like (-) part is dominantly controlled by the width 
7 - relative to the pole mass - which determines the fraction of Tqq with negative 
invariant mass squared (P 2 < 0) relative to the time-like part Tq Q . 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 1 2 3 4 5 



p [GeV/c] p [GeV/c] 

Fig. 4. The integrand I(u,p) (16) for T = 1.02T C (l.h.s.) and T = 2T C (r.h.s.) in terms of 
contour lines. The straight (blue) line (u = p) separates the lime-like from the space-like 
sector. Note that for a convergence of the energy density integral the upper limits for oj 
and p have to be increased by roughly an order of magnitude compared to the area shown 
in the figure. 

The explicit results for the quasiparticle energy densities T j) and T q are displayed 
in Fig. 5 by the dashed blue and dot-dashed red lines (multiplied by (Tc/T) 4 ), 
respectively. As in case of N + and N~ the space-like energy density Tqq" is seen 
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Fig. 5. Upper part: The time-like energy density TJ^ (dashed blue line), the space-like energy 
density Tqq" (dot-dashed red line) and the total energy density Too = Tqq + (thick solid 
green line) as a function of T/T c . The thin black line displays the energy density e(T/T c ) 
from (9); it practically coincides with Too within the linewidth and is hardly visible. All 
densities are multiplied by the dimensionless factor (T c /T) 4 in order to divide out the 
leading temperature dependence. Lower part: Same as the upper part in order to enhance 
the resolution close to T c . 

to be larger than the time-like part T ^ for all temperatures above 1.05 T c . Since 
the time-like part corresponds to the independent quasiparticle energy density 
within the lightcone, the space-like part Tqq" can be interpreted as an interaction 
density V if the quasiparticle energy T 00 matches the total energy density e(T) (9) 
as determined from the thermodynamical relations (8) and (9). In fact, the DQPM 
yields an energy density T o - adding up the space-like and time-like parts - that 
almost coincides with e(T) from (9) as seen in Fig. 5 where both quantities (multiplied 
by (Tc/T) 4 ) are displayed in terms of the thin black and thick solid green lines, 
respectively; actually both results practically coincide within the linewidth for T > 
2T C . An explicit representation of their numerical ratio gives unity within 2% for 
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T > 2T C ; the remaining differences can be attributed to temperature derivatives 
~ d/dT(ln('~f / E)) etc. in order to achieve thermodynamic consistency but this is 
not the primary issue here and will be discussed in a forthcoming study [42]. The 
deviations are more clearly visible close to T c (lower part of Fig. 2) where the variation 
of the width and mass are most pronounced. However, for all practical purposes one 
may consider T 00 (T) ps e(T) and separate the kinetic energy density TjJ from the 
potential energy density T q as a function of T or - in equilibrium - as a function of 
the scalar gluon density N s or iV + , respectively. 



3 Dynamics of time-like quasiparticles 

Since in transport dynamical approaches there are no thermodynamical Lagrange 
parameters like the inverse temperature (3 = T _1 or the quark chemical potential 
fi q , which have to be introduced in thermodynamics in order to specify the average 
values of conserved quantities (or currents in the relativistic sense), derivatives of 
physical quantities with respect to the scalar density p s = N s (or time-like gluon 
density p g = N + ) are considered in the following (cf. Ref. [43]). As mentioned above 
one may relate derivatives in thermodynamic equilibrium via, 



if the volume and pressure are kept constant. For example, a numerical evaluation 



with b= 5, ai = 1.5/m~ 3 and a 2 = 104/m~ 3 , which follows closely the quadratic 
scaling in T/T c as expected in the Stefan Boltzmann limit. The additional exponential 
term in (18) provides a sizeable correction close to T c . The approximation (18) may 
be exploited for convenient conversions between p s and T/T c in the pure gluon case 
but will not be explicitly used in the following. 

The independent quasiparticle energy density T K := T ~q and potential energy density 
V := Tqq now may be expressed as functions of p s (or p g ) instead of the temperature 
T. The interaction energy density then might be considered as a scalar energy density 
which - as in the nonlinear a-model for baryonic matter [44] - is a nonlinear function 
of the scalar density p s . As in case of nuclear matter problems the scalar density 
p s does not correspond to a conserved quantity when integrating over space; it only 
specifies the interaction density parametrically, i.e. V(p s ). Alternatively one might 
separate V into parts with different Lorentz structure, e.g. scalar and vector parts 



d d dp s 
dT ~ dp~ s ~dT' 



(17) 



of dp s /d(T/T c ) gives 




(18) 
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as in case of nuclear matter problems [44], but this requires additional information 
that cannot be deduced from the DQPM alone. 
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Fig. 6. Upper part: The quasiparticle energy per degree of freedom Tk/N + (dashed blue 
line) and the space-like potential energy per degree of freedom V/N + (dot-dashed red line) 
as a function of T/T c . All energies are multiplied by the dimensionless factor (T c /T). Lower 
part: Same as the upper part in order to enhance the resolution close to T c . 



It is instructive to show the 'quasiparticle' and potential energy per degree of free- 
dom Tk/N + and V/N + as a function of e.g. N + , N s or T/T c . As one might have 
anticipated the kinetic energy per effective degree of freedom is smaller than the 
respective potential energy for T/T c > 1.05 as seen from Fig. 6 where both quan- 
tities are displayed as a function of T/T c in terms of the dashed and dot-dashed 
line, respectively. It is seen that the potential energy per degree of freedom steeply 
rises in the vicinity of T c whereas the independent quasiparticle energy rises almost 
linearly with T. Consequently rapid changes in the density - as in the expansion 
of the fireball in ultrarelativistic nucleus-nucleus collisions - are accompanied by a 
dramatic change in the potential energy density and thus to a violent acceleration of 
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the quasi-particles. It is speculated here that the large collective flow of practically 
all hadrons seen at RHIC [6] might be attributed to the early strong partonic forces 
expected from the DQPM. 




OLd , i , i 

1 10 100 
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Fig. 7. The mean-field potential U(N + ) = U(p g ) as a function of the time-like gluon density 
iV + = p g in comparison to the fit (19) (solid blue line). The densities N + = 1, 1.4, 5, 10, 
50, 100 fm -3 correspond to scaled temperatures of T/T c ~ 1.025, 1.045, 1.25, 1.5, 2.58, 
3.25, respectively (cf. Fig. 3). 

In order to obtain some idea about the mean-field potential U s (p s ) (or U(p g ) in the 
rest frame) one can consider the derivative dV/p s = U s (p s ) or dV/N + = U(N + ) = 
U(p g ). The latter is displayed in Fig. 7 as a function of N + = p g and shows a distinct 
minimum at p g ~ 1.4 fm~ 3 which corresponds to a temperature T pa 1.045T C . The 
actual numerical results can be fitted by the expression, 

U( Pg ) = 4^-39 e -^/°- 31 + 2.93 p° g 21 + 0.55 pf 6 [GeV] , (19) 

dPg 



where p g is given in fm -3 and the actual numbers in front carry a dimension in order 
to match to the proper units of GeV for the mean-field U. By analytical integration 
of (19) one obtains a suitable approximation to V(p g ). The approximation (19) works 
sufficiently well as can be seen from Fig. 7 - showing a comparison of the numerical 
derivative dV/dN + with the fit (19) in the interval 0.7 fm~ 3 < N + < 300 fm~ 3 
- such that one may even proceed with further analytical calculations. Note that 
a conversion between the time-like quasiparticle density N + = p g and the scalar 
density p s is easily available numerically (cf. lower part of Fig. 3) such that derivatives 
with respect to p s are at hand, too; the latter actually enter the explicit transport 
calculations [45] while derivatives with respect to p g in the rest frame of the system 
are more suitable for physical interpretation and will be used below. 
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Some information on the properties of the effective gluon-gluon interaction v gg may 
be extracted from the second derivative of V with respect to p g , i.e. 

VgM := "TT ~ " 125 - 8 e ~ P9/ °' 31 + 0.615/pJ 79 + 0.2/p°- 64 [GeVfm 3 ], (20) 
dp g 



where the numbers in front have again a dimension to match the units of GeV 
fm 3 . The effective gluon-gluon interaction v gg (20) is strongly attractive at low den- 
sity 0.003 fm~ 3 < p g and changes sign at p g ~ 1.4 fm~ 3 to become repulsive at 
higher densities. Note that the change of quasiparticle momenta (apart from col- 
lisions) will be essentially driven by the (negative) space-derivatives —VU(x) = 
—dU(pg) I 'dp g V p g (x) (or alternatively by —dU s (p s )/dp s Vp s (x)). This implies that 
the gluonic quasiparticles (at low gluon density) will bind with decreasing density, 
i.e. form 'glueballs' dynamically close to the phase boundary and repell each other 
for p g > 1.4 fm~ 3 . Note that color neutrality is imposed by color-current conser- 
vation and only acts as a boundary condition for the quantum numbers of the 
bound/resonant states in color space. 

This situation is somehow reminiscent of the nuclear matter problem [44] where a 
change in sign of the 2nd derivative of the potential energy density of nuclear matter 
at low density indicates the onset of clustering of nucleons, i.e. to deuterons, tritons, 
a-particles etc., which form the states of the many-body system at low nucleon 
densities (and not a low density nucleon gas). This is easy to follow up for the 
simplified nonrelativistic energy density functional e N for nuclear matter, 

e N - Apf + \p\ + fpf, (21) 



where the first term gives the kinetic energy density and the second and third term 
correspond to attractive and repulsive interaction densities. For A ps 0.073GeV fm 2 , 
B w —1.3 GeV fm 3 and C ps 1.78 GeV fm 4 a suitable energy density for nuclear 
matter is achieved; it gives a minimum in the energy per nucleon E/A = e^/p^ ~ 
—0.016 GeV for nuclear saturation density p° N ~ 0.168 fm~ 3 . The mean-field potential 
Un = Bpw+Cp 4 ^ 3 has a minimum close to p° N such that the effective nucleon-nucleon 

1 /3 

interaction strength vnn = B + 4/3Cp^ changes from attraction to repulsion at 
this density. Note that in the gluonic case the minimum in the mean-field potential 
U (19) occurs at roughly 8 times p% and the strength of the gluonic interaction is 
higher by more than 2 orders of magnitude! 

The confining nature of the effective gluon-gluon interaction v gg (20) becomes ap- 
parent in the limit p g — > 0, where the huge negative exponential term dominates for 
p g > 0.003 fm~ 3 ; for even smaller densities the singular repulsive terms take over. 
Note, however, that the functional extrapolation of the fit (19) to vanishing gluon 
density p g has to be considered with care and it should only be concluded that the 
interaction strength becomes 'very large'. On the other hand the limit p g — > is only 
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academical because the condensation/fusion dynamically occurs for p g ps 1 fm 3 . 

A straight forward way to model the gluon condensation or clustering to confined 
glueballs dynamically (close to the phase transition) is to adopt a screened Coulomb- 
like potential v c (r, A) with the strength / d 3 r v c (r, A) fixed by v gg (p g ) from (20) and 
the screening length A from 1QCD studies. For the 'dilute gluon regime' (p g < 1.4 
fm -3 ), where two-body interactions should dominate, one may solve a Schrodinger 
(or Klein-Gordon) equation for the bound and/or resonant states. This task is not 
addressed further in the present study since for the actual applications (as in the 
Parton-Hadron-String-Dynamics (PHSD) approach [45]) dynamical quark and anti- 
quarks have to be included. The latter degrees of freedom do not change the general 
picture very much for higher temperatures T > 2T C but the actual numbers are dif- 
ferent close to T c since the quarks and antiquarks here dominate over the gluons due 
to their lower mass. The reader is referred to an upcoming study in Ref. [42]. 

Some comments on expanding gluonic systems in equilibrium appear in place, i.e. 
for processes where the total volume V and pressure P play an additional role. For 
orientation we show the entropy per time-like particle s/N + in Fig. 8 as a function of 
A^ + (upper) and T/T c (lower part) which drops close to the phase boundary since the 
quasiparticles become weakly interacting (cf. Fig. 6). Note that this is essentially due 
to the low density and not due to the interaction strength (20); a decrease of the width 
7 (as encoded in (5)) implies a decrease in the interaction rate! An expansion process 
with conserved total entropy S = sV leads to a change in the total gluon number 
N + V since s/N + changes with density (or temperature) (Fig. 8). The same holds 
for an expansion process with constant total energy eV since also e/N + is varying 
with density (or temperature). Other scenarios involving e.g. S = P/T also involve 
a change of the gluon number N + V during the cooling process such that reactions 
like gg <-> g, ggg <-> gg etc. are necessary ingredients of any transport theoretical 
approximation. We do not further investigate different expansion scenarios here since 
the reactions g <-> qq, i.e. the gluon splitting to a quark and antiquark as well as 
the backward fusion process, are found to play a dominant role in the vicinity of the 
phase transition as well as for higher temperatures [42,45]. 



4 Conclusions and discussion 

The present study has provided a novel interpretation of the dynamical quasiparticle 
model (DQPM) by separating time-like and space-like quantities for 'particle densi- 
ties', energy densities, entropy densities ect. that also paves the way for an off-shell 
transport approach [45]. The entropy density s in (7) is found to be dominated by the 
on-shell quasiparticle contribution (first line in (7)) (cf. [31]) while the space-like part 
of the off-shell contribution (second line in (7)) gives only a small (but important) 
enhancement (cf. Fig. 2). However, in case of the 'gluon density' N = A^ + + N~ and 
the gluon energy density T 0Q = T^ + T^ the situation is opposite: here the space-like 
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Fig. 8. The entropy per degree of freedom s/N + as a function of N + (upper part) or T/T c 
(lower part). 

parts (A^",T q) dominate over the time-like parts (N + ,Tq ) except close to T c where 
the independent quasiparticle limit is approximately regained. The latter limit is a 
direct consequence of the infrared enhancement of the coupling (4) close to T c (in 
line with the 1QCD studies in Ref. [36] ) and a decrease of the width 7 (5) when 
approaching T c from above. 

Since only the time-like part iV + can be propagated within the lightcone the space- 
like part N~ has to be attributed to t-channel exchange gluons in scattering processes 
that contribute also to the space-like energy density T q. The latter quantity may 
be regarded as potential energy density V. This, in fact, is legitimate since the 
quasiparticle energy density T 00 very well matches the energy density (9) obtained 
from the thermodynamical relations. Only small deviations close to T c indicate that 
the DQPM in its straightforward application is not thermodynamically consistent. 
However, by accounting for 'rearrangement terms' in the energy density - as known 
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from the nuclear many-body problem [46] - full thermodynamical consistency may 
be regained [42]. 

It is instructive to compare the present DQPM to other recent models. In the 
PNJL 2 model [47] the gluonic pressure is build up by a constant effective potential 
£/($,$*; T) which controls the thermodynamics of the Polyakov loop $. It is ex- 
panded in powers of $$* with temperature dependent coefficients in order to match 
1QCD thermodynamics. Thus in the PNJL there are no time-like gluons; the effec- 
tive potential U(Q, <3>*;T) stands for a static gluonic pressure that couples to the 
quark/ ant iquark degrees of freedom. The latter are treated in mean-field approxi- 
mation, i.e. without dynamical width, whereas the DQPM incorporates a sizeable 
width 7. 

Another approach to model 1QCD thermodynamics has been suggested in Ref. [43] 
and is based on an effective Lagrangian which is nonlinear in the effective quark and 
gluon fields. In this way the authors avoid a parametrization of the interaction density 
in terms of Lagrange parameters (T, fi) and achieve thermodynamical consistency. 
The latter approach is closer in spirit to the actual interpretation of the DQPM and 
may be well suited for an on-shell transport theoretical formulation. The on-shell 
restriction here comes about since effective Lagrangian approaches should only be 
evaluated in the mean-field limit which implies vanishing scattering width for the 
quasiparticles. This is sufficient to describe systems is thermodynamical equilibrium, 
where forward and backward interaction rates are the same, but might not provide 
the proper dynamics out-of-equilibrium. 

Some note of caution with respect to the present DQPM appears appropriate: the pa- 
rameters in the effective coupling (4) and the width (5) have been fixed in the DQPM 
by the entropy density (7) to 1QCD results assuming the form (2) for the spectral 
function p{uj). Alternative assumptions for p{uj) will lead to slightly different results 
for the time-like density, energy densities etc. but not to a qualitatively different 
picture. Independent quantities from 1QCD should allow to put further constraints 
on the more precise form of p(u) such as calculations for transport coefficients [16]; 
unfortunately such 1QCD studies are only at the beginning. A more important issue 
is presently to extend the DQPM to incorporate dynamical quark and antiquark 
degrees of freedom (as in [31]) in order to catch the physics of gluon splitting and 
quark- ant iquark fusion (g^q + q, g + g^q + q + g) reactions [42,45]. 

Coming back to the questions raised in the Introduction concerning i) the appropri- 
ate description of QCD thermodynamics within the DQPM and ii) the possibility to 
develop a consistent off-shell partonic transport approach as well as iii) the perspec- 
tives for a dynamical description of the transition from partonic to hadronic degrees 
of freedom, we are now in the position to state: most likely 'Yes'. 



2 Polyakov-loop-extended Nambu Jona-Lasinio 
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